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Abstract 

Eigenvalues and eigenfunctions of the QCD Dirac operator are studied for gauge field 
configurations given by a liquid of instantons. We find that for energy differences 5E below 
an energy scale E c the eigenvalue correlations are given by Random Matrix Theories with 
the chiral symmetries of the QCD partition function. For eigenvalues near zero this energy 
scale shows a weak volume dependence that is not consistent with E c ~ 1/L 2 which might 
be expected from the pion Compton wavelength and from the behavior of the Thouless 
energy in mesoscopic systems. However, the numerical value of E c for our largest volumes 
is in rough agreement with estimates from the pion Compton wavelength. A scaling 
behaviour consistent with E c ~ 1/L 2 is found in the bulk of the spectrum. For SE > E c 
the number variance shows a linear dependence with a slope which is larger than the 
nonzero multifractality index of the wave functions. Finally, the average spectral density 
and the scalar susceptibilities are discussed in the context of quenched chiral perturbation 
theory. We argue that a nonzero value of the disconnected scalar susceptibility requires 
a linear dependence of the number variance on 5E. 



1 Introduction 



Random matrix theories have been applied extensively to disordered mesoscopic systems 
(for recent reviews we refer to Jl|, |3|). In particular, correlations of energy eigenvalues 
have been studied in great detail. One important energy scale that enters in these studies 
is the Thouless energy, E c , which is essentially the inverse diffusion time of an electron 
through the sample |4j] . If the diffusion constant is given by D and the linear dimension 
of the sample is equal to L the Thouless energy is given by 
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What has been found is that eigenvalues closer than E c are correlated according to the 
invariant random matrix ensembles. Because such correlations are characteristic for classi- 
cally ergodic and chaotic quantum system this regime is also known as the ergodic energy 
domain. A second energy scale is given by h/r e , where r e is the elastic collision time. 
The energy regime with 5E > h/r e is known as the ballistic regime whereas the domain 
E c < 5E < h/r e is known under the names of diffusive regime or Altshuler-Shklovskii 
regime. The different energy regimes have received a great deal of attention in recent 
literature ]5|, |5], 0, |S], [II], |TT[ (for further references we refer to |TJ, [| [3|]). A semiclassical 
interpretation in terms of the return probability has been given in jl2| . 

A convenient way to measure the pair correlations of eigenvalues is by means of the 
number variance S 2 (n). This statistic is defined as the variance of the number of eigenval- 
ues in an interval containing n eigenvalues on average. In the ergodic regime the number 
variance is given by the invariant random matrix ensembles £2(77.) ~ (2//37T 2 ) log(n). The 
behavior of the number variance in the diffusive is less well-established. One proposal is 
@] that S 2 (n) ~ n d l 2 (where d is the number of Euclidean dimensions) if the coupling is 
below the critical value for a localization transition. However, for critical values of the 
coupling constant, for which the localization length scales with the size of the sample, a 
linear dependence on n is predicted, i.e., ^(n) ~ \n. In this regime, the slope has been 



related to the multifractality index of the wave functions |T3|, [14 . 

In this letter we investigate to what extent such scenarios are realized in QCD. In 
QCD the order parameter of the chiral phase transition is related to the spectral density, 



p(A), of the Dirac operator by means of the Banks-Casher formula [15|, X = Iim7rp(0)/V 
(the space-time volume is denoted by V = L 4 ). For broken chiral symmetry, E ^ 0, this 
results in eigenvalues near zero virtuality that are spaced as ~ 1/V". It is therefore natural 
to introduce the microscopic limit V — > oo such that u = WE is kept constant. The 
microscopic spectral density is then defined by ]|16] , 
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where the spectral density is defined by 
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and the brackets denote ensemble averaging. There is ample evidence from lattice QCD 
17, 18, 19] and instanton liquid simulations |20 that ps{u) and other correlators on 



the scale of individual level spacings |2T| are given by chiral Random Matrix Theory, 
i.e. random matrix theories with the chiral symmetries of the QCD partition function. 



However, at scales beyond a few eigenvalue spacings in both instanton simulations |20| and 
lattice QCD simulations |17|, IS] the Dirac eigenvalues near zero show stronger fluctuations 
than those in the chiral random matrix theories. This indicates the presence of an energy 
scale in QCD which may be identified as the Thouless energy. 

In the regime where the pion loops can be ignored the correlations of the eigenvalues 
near zero are given by the chiral random matrix ensembles E3J . They can be mapped 
onto an effective partition function given by the mass term generated by the spontaneous 
breaking of chiral symmetry ]TE| . From the work of Gasser, Leutwyler and Smilga [23, |23fl, 
we expect that the boundary of the ergodic regime is given by a mass scale m c where the 
Compton wavelength of the associated Goldstone boson with mass m 2 = m c K (according 
to the PCAC relation K = E/F 2 , where F n is the pion decay constant) is of the order of 



the linear dimension of the box. This condition can be written as |25|, [26 
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It is therefore tempting to interpret 1/K as the diffusion constant. Because K is large 
on a hadronic scale (K 1660 MeV) the diffusion constant is relatively small, and we 
expect a Thouless energy that is small as well. We emphasize that m c is a valence quark 



mass that is not present in the fermion determinant. In other words, the spectral density 
does not depend on m c . For example, the valence quark mass of the chiral condensate is 
related to the spectral density by Wi 



%)4f^. (5) 
V Jo X z + m z 

Therefore, the spectral density and, similarly, higher order correlation functions, can be 
obtained from the valence quark mass dependence of the effective partition function. It 
can be obtained with the help of quenched chiral perturbation theory |28|] or by means 
of the replica trick in which we start from n/ valence flavors and take the limit n/ — > 
at the end of the calculation. In the thermodynamic limit, the mass scale (^) is much 
less than the scale of massive quarks which leads to a decoupling of the valence quark 
mass dependent part of the partition function. The Dirac spectrum and its correlations 
are thus given by the valence quark mass dependence of the effective partition function 
which below the scale @ is given by chiral random matrix theory. For massless quarks 
the argument is simpler. Now the Compton wavelength of the physical pions is always 
much larger than the size of the box. The crossover point is thus given by the value of the 
valence quark mass where the corresponding Goldstone boson (in the sense of quenched 
or partially quenched chiral perturbation theory) is equal to the size of the box. In this 
case, however, the microscopic spectral correlations depend on the number of massless 
flavors. Alternatively, one can study the boundaries of the domain where the pion loops 
can be neglected via relations between microscopic spectral correlators and finite volume 



partition functions as given in [29 



In the bulk of the spectrum, chiral symmetry is broken explicitly by the valence quark 
mass, and the 'Goldstone bosons' associated with the valence quark masses are no longer in 
the low energy spectrum of QCD. However, in this case, spontaneous symmetry breaking 
in the generating function for the correlation functions results in the usual nonlinear a— 
model for disordered mesoscopic systems [50]. The critical energy scale is then given 



by the Thouless energy m c ~ 1/K'L 2 , where K' = irp(E)/VF% and p(E) is the spectral 
density at E. Since the group manifold is different in this case |H[ the numerical constant 
in this relation differs by a factor of the order unity from the constant in ([|). 

Using that according to the Banks Casher relation the eigenvalue spacing is given 



by A = 7r/EV, the condition (f|) can be rewritten in dimensionless form as 

A 7T 

Here, g c plays the role of the dimensionless conductivity. At this point we wish to mention 
that a relation between transport properties of mesoscopic systems and the pion decay 
constant was established recently in [32]. In lattice QCD for an Na A lattice this relation 
reads g c = F^a 2 \^N /n. On a 16 4 lattice with a lattice spacing of 0.1 fm this results in a 
dimensionless Thouless energy of a couple of lattice spacings. 

In order to investigate the presence of a Thouless energy in QCD we perform our 
calculations for gauge field configurations given by a liquid of instantons. For such con- 
figurations we study correlations of the eigenvalues of the corresponding Dirac operator. 
In particular, we investigate the number variance, and establish a crossover between the 
ergodic and the diffusive domain. The linear dependence will be compared with the mul- 
tifractality index of the wave functions. We will argue that results for chiral perturbation 
theory apply to the diffusive regime. The above observables are studied for a variety of 
system sizes and as a function of the number of flavors. 



2 Instanton liquid model 

In this model the gauge field configurations are given by a superposition of instantons. We 
use the streamline Ansatz which for a dilute system reduces to the sum Ansatz defined as 
a simple superposition of instanton profiles = J2i Az>- Each instanton is described by 
4N C (the number of colors is denoted by N c ) collective coordinates. The Euclidean QCD 
partition function is then approximated by 

Z- mst = J DO. det^( 7 • D + m)e" SYM , (7) 

where the integral is over the collective coordinates of the instantons and the covariant 
derivative is defined by = + iA^. The Yang- Mills action of the gauge field configu- 
rations is denoted by SVm- We will only consider configurations with an equal number of 
instantons and anti-instantons. The fermion determinant is evaluated in the space of the 
fermionic zero modes of the instantons. 
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This partition function obeys the flavor and chiral symmetries of the QCD partition 
function. It reproduces results obtained in chiral perturbation theory as well as many other 
details of the low energy hadronic phenomenology |33| , [34], [35], |36| . In our calculations we 
use the standard instanton liquid parameters with instanton density N/V = 1 (in units 
of fmT 1 ). For more details on the instanton liquid partition function we refer to a recent 
review by Schafer and Shuryak [[37j . 

3 Random matrix model 

The chiral random matrix ensembles for Nf massless quarks in the sector of topological 



charge v are defined by the partition function |16, |22| 



= J DWd*»' ( ^ ™ J (8 ) 

where W is a n x (n+v) matrix. As is the case in QCD, we assume that v does not exceed 
\/2n. The parameter 2n is identified as the dimensionless volume of space time. It can 
be thought of as the total number of instantons in the partition function (0). The matrix 
elements of W are either real {(3=1, chiral Orthogonal Ensemble (chOE)), complex 
{(3 = 2, chiral Unitary Ensemble (chUE)), or quaternion real {(3 = 4, chiral Symplectic 
Ensemble (chSE)). The simplest ensemble is the Gaussian case with V{x) = Y?x (also 
known as the Laguerre ensemble). In that case one can easily show that the microscopic 
spectral density (0) is given by |22| 



Ps{u) = ~{J 2 a {u) - J a+1 {u)J a ^{u)). (9) 
where a = Nf + \u\ (the result for a = was obtained in |38||). It is also straightforward 



to derive the microscopic correlation functions |38], |18|| . In the bulk of the spectrum, the 
eigenvalue correlations of the chiral ensembles are given by the invariant random matrix 
ensembles. 



It was shown by Akemann et al. |55 that, for (3 = 2, the microscopic spectral density 
and the microscopic spectral correlators do not depend on the potential V{x) and are 
given by the results PSI for the Laguerre ensemble. Their work greatly extends an earlier 



result by Brezin and Zee ||D[ who studied the potential V{x) = x + ax . The microscopic 



spectral density is also invariant with respect to the addition of an arbitrary constant 



matrix to the off-diagonal blocks in the determinant |4l], [42| . The application of chRMT 
has been put on a firm foundation by these and other universality proofs []43|| . Whether 
or not QCD is in this universality class is a dynamical question that only can be proven 
by explicit numerical simulations. 

4 Numerical Results 

In this section we simulate the instanton liquid model introduced in section 2. We evaluate 
the partition function (|7]) by means of a Metropolis algorithm. Typically, we perform on 
the order of 10,000 sweeps for each set of parameters. For the equilibrated configurations 
the eigenvalues and eigenvectors of the Dirac operator are calculated by means of standard 
diagonalization procedures. Most of our calculations have been performed in the quenched 
approximation (Nf = 0). The motivation for this choice is twofold. First, calculations 
at Nf ytz o require the evaluation of a determinant for each update of the integration 
variables which is extremely costly and restricts us to a total number of instantons not 
larger than 128. For Nf = 0, we are able to generate ensembles up to 512 instantons 
which turns out to be important in the study of the thermodynamic limit. Second, our 
results are compared with ideas from disordered mesoscopic systems where no fermion 
determinant is present. 

In order to separate the fluctuations of the eigenvalues from the average spectral den- 
sity we unfold the spectrum, i.e. we rescale the spectrum in units of the local average level 
spacing. Specifically, the unfolded spectrum A™ f is obtained from the original spectrum 
Xk according to 

A£ nf = f Xk p(s)ds, (10) 
Jo 

where p(s) is the smoothened average spectral density. The number variance, the nearest 
neighbor spacing distribution and the microscopic spectral density are all calculated from 
the unfolded spectrum. 



In Figs. 1 and 2 we show the number variance, ^(n) versus n for Nf = and various 
total numbers of instantons as given in the legend of the figure. The instanton density 
N/V = 1. The random matrix result for E 2 (n) for the chGUE |18| is depicted by the 
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Fig. 1. The number variance ^(n) versus n in the quenched approximation 
for an interval starting at A = 0. The total number of instantons is indicated 
in the legend of the figure. 

solid curve. In both figures, the right figure (with results for only two different volumes) 
is a blown-up version of the left figure. In Fig. 1 the number variance is calculated for the 
interval that contains on average the n smallest positive eigenvalues. Fig. 2 represents the 
number variance in the bulk of the spectrum obtained from an interval that is symmetric 
about the average unfolded positive eigenvalue. In both cases we observe a transition 
point n c below which the number variance is given by random matrix theory. In Fig. 1 
the value of the crossover point, n c ^ 2, depends only weakly on the total number of 
instantons (or the volume). This is not in agreement with the theoretical expectation 
([]) that n c m F^y/V /ir in four dimensions. For standard instanton liquid parameters 
the numerical value of n c for N instantons is given by n c « O.OSv^V which is on the 
order of the results found in Fig. 1. In the bulk of the spectrum (Fig. 2) the value of 
n c is consistent with a \/V scaling but the numerical constant appears to be larger than 
the above estimate (even if we take into account that the spectral density in the bulk 
is considerably less than near A = 0). This result is in agreement with the finding that 
correlations of lattice QCD eigenvalues for 12 4 and 16 4 lattices are given by RMT up to 
distances of 100 and 300 level spacings, respectively |21, 44|. Notice that the condition 
that the pion Compton wavelength is equal to the size of the box (leading to (f|)) provides 
us only with an order of magnitude estimate for the transition point. 
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Fig. 2. The number variance ^(n) versus n in the bulk of the spectrum in 
the quenched approximation and a total number of instantons as indicated in 
the legend of the figure. 

Beyond the crossover point the number variance shows a linear behavior with a slope 
X ~ 0.08 for eigenvalues near zero and x ~ 0.04 in the bulk of the spectrum. The 
downward trend of the curves for larger values of n is a finite size effect. For example, for 
a finite number of uncorrelated eigenvalues S 2 (n) does not behave linearly, but rather as 
S 2 (n) = n — 2n 2 /N (notice that we have only N/2 positive eigenvalues). This finite size 
effect prevents us from saying more about the ballistic regime, an energy scale of roughly 
the inverse distance between instantons. 

In the ergodic regime we expect that eigenvalue correlations are given by the invariant 
random matrix ensembles. This is shown in Fig. 3 where we plot the nearest neighbor 
spacing distribution P(S) versus S (left) and the microscopic spectral density ps( z ) versus 
z (right) (in terms of the variable z = u/ir in (Q)). The full circles and the open squares 
represent results for 512 and 128 instantons, respectively. Results obtained for the (chiral) 
random matrix ensembles are depicted by the full line. Clearly, below the Thouless energy 
we are in close agreement with RMT. However, beyond the Thouless energy of about two 
level spacings the oscillatory structure in the microscopic spectral density is no longer 
reproduced. This implies that beyond this point the QCD Dirac eigenvalues fluctuate 
more than the chGUE eigenvalues which is consistent with the number variance shown 
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Fig. 3. The nearest neighbor spacing distribution P(S) versus S (left) and 
the microscopic spectral density ps(z) versus z (right). The full circles and 
the open squares represent results for 512 and 128 instantons, respectively. 
Results obtained from the (chiral) random matrix ensembles are depicted by 
the full line. 

in Fig. 1. A similar increase in eigenvalue fluctuations at a scale of a few eigenvalue 



spacings has been observed in lattice QCD spectra [17]] indicating that also in that case 
the Thouless energy for the region around A = is less than predicted by the scaling 
behavior of disordered mesoscopic systems. 

The inverse participation ratio which is a measure for the number of significant com- 
ponents in the wave function is defined as 

/ 2 (A) = (£hMA)l 4 >, (11) 

k 

The V'fc(A) are the normalized eigenf unctions of the Dirac operator in the space of fermionic 
zero modes corresponding to eigenvalue A. The number of components of the wave func- 
tion will be denoted by N. For the chGUE we find I2 = 2/(N+2) whereas for uncorrelated 
eigenvalues (the chiral Poisson ensemble) we find I2 = 1/2. Results for Nl2(X) as a func- 
tion of A are depicted in Fig. 4 (left). The general impression is that the wave-functions 
are extended with a participation ratio that is not too different from the random ma- 
trix result (full line). The eigenf unctions corresponding to small and large eigenvalues 
appear to be somewhat more localized. A more definitive result for the character of the 
wavefunctions follows from the scaling behavior of /2(A) with the volume. 
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Fig. 4. The inverse participation ratio times N versus the corresponding Dirac 
eigenvalues for Nf = (left) and a total number of instantons equal to 64, 
128, 256 and 512, and the scaling behavior versus N (right) for the energy 
intervals [0.1,0.2] and [0.7,0.8]. 

A double logarithmic plot of Nl2(X) versus N is shown in Fig. 4 (right). Results are given 
both for the energy intervals [0.1,0.2] (open circles) and [0.7,0.8] (full circles). The first 
region corresponds to a part of the spectrum where the number variance shows a linear 
behavior for the volumes shown in Fig. 1, and the second region corresponds to the bulk 
of the spectrum. 

The multifractality index rj of the wave functions is denned by [13 



yv/d- 



(12) 



According to an argument given in []jl| the value of rj = 2\d (where x 1S t ne slope of the 
linear piece in the number variance) in the critical domain. From the volume dependence 
of the inverse participation ratio shown in Fig. 4 (right) we find a value for r\jd of about 
0.02 and 0.04 for the intervals [0.1,0.2] and [0.7,0.8], respectively. These values are well 
below the theoretical result of 1%. Apparently, the ensemble of instantons is not in the 
critical region for Nf = 0. The localization properties of eigenfunctions have also been 
studied for the Wilson lattice QCD Dirac operator ||45|| . In that case it was found that 
the eigenstates are localized. We have no explanation for this discrepancy. 

The flavor dependence of the number variance of eigenvalues near zero is shown in Fig. 
5. Results are given for 128 instantons and massless flavors. Again the Thouless energy 
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is about two level spacings (in order to avoid a cluttered figure we have not plotted 
the random matrix results for the number variance). We observe an enhancement of 
eigenvalue fluctuations for 3 flavors (chiral symmetry is already restored for Nf = 4 f37| ) 
and conclude that the critical number of flavors is approximately equal to three. However, 
our volumes (of up to 128 instantons) were too small to reliably extract the multifractality 
index. For completeness we mention that a similar enhancement of fluctuations for Nf = 3 
also takes place in the bulk of the spectrum. 




n 



Fig. 5. The flavor dependence of the number variance ^(ji) versus n calcu- 
lated for an ensemble of 128 instantons. 

5 Spectral density and chiral susceptibility 



In this section we consider the average spectral density, the valence quark mass depen- 
dence of the chiral condensate and the connected and disconnected scalar susceptibilities. 
In order to compare our results with quenched chiral perturbation theory |28[ we put 
particular emphasis on the extrapolation to the thermodynamic limit. These quantities 
were calculated elsewhere for smaller volumes (64 instantons) p6l [46]. 
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The spectral density of the Dirac eigenvalues for a liquid of instantons was studied 
before for Nf = with the help of mean field techniques |34|, [33|. For large values of iV c 



the overlap matrix elements are independent resulting in a semi-circular distribution |34 
For small values of N c the matrix elements are correlated and the natural distribution 
is a Gaussian [33]. On the other hand, from the asymptotic dependence of the matrix 



elements of the Dirac operator it was concluded that ||71 the spectral density diverges as 
1/y/X. Another estimate comes from quenched chiral perturbation theory [p8 |. For A — > 



one obtains a linear dependence |48[ for Nf ^ 0, and a logarithmic divergence j49| for 
Nf = 0. For Nf the slope of the Dirac spectral density for A — > is consistent with 
the analytical result [ 5Q(| . Here we only consider the quenched case. 



In order to study the thermodynamic limit of the spectral density we have unfolded 
the calculated spectrum with the known result for the microscopic spectral density. The 
resulting normalized average spectral density for Nf = is shown in Fig. 6. Inspection of 
the normalized eigenvalue density in an interval [0.0, 0.05] shows that it grows logarithmi- 
cally with N. This suggests that p(A) diverges for A — > in the thermodynamic limit. In 
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order to understand the nature of this divergence we have compared results obtained for 
the instanton liquid with the analytical results for the a chiral version of the a-ensenbles 
|5lH (dotted curve). These are random matrix ensembles with potential V(x) = x a ^ 2 in 
(§). The spectral density of the chiral a-ensembles diverges logarithmically for A — > at 
a — 1. The value of p(0) is finite for a > 1 . For iV = 64, 128, 256 and 512 we find 
a = 1.217, 1.132, 1.073, 1.039 by fitting to the energy interval [0.05,0.5], respectively. In 
the thermodynamic limit this extrapolates to a value of a = 1.00 ±0.02, resulting in a log- 
arithmic divergence of the spectral density. We also compare the average spectral density 
to a simple power law behavior a — bX 13 . We find that the spectra are equally well fitted 
by such dependence (see full curve in Fig. 6). By fitting to the interval [0.05, 0.5] we find 
(3 = 0.28, 0.25, 0.13, 0.09 for N = 64, 128, 256, and 512, respectively. This extrapolates 
to a value of (3 = 0.03 ± 0.04 and is also consistent with a logarithmic singularity for 
A — > 0. We conclude that our results are consistent with an average spectral density that 
diverges logarithmically for A — > in the limit V — > oo. However, larger volumes need 
to be studied in order to arrive at definite conclusions. We wish to emphasize that our 
results disagree with estimates obtained by mean field methods. 

The valence quark mass dependence of the chiral condensate defined in (|5|) is shown 
in Fig. 7. Results are given for Nf = and different values of the total number of 
instantons. The downward trend for small masses is a finite size effect that can be 
understood analytically in terms of the microscopic spectral density [^5] . The large volume 
dependence of the "plateau" prevents us from comparing our results with quenched chiral 



perturbation theory |28[ which predicts a logarithmic dependence on the valence quark 
mass. For the same reason, the connected scalar susceptibility, <9 m £(m), shows an even 
larger volume dependence. 

Finally, we consider the disconnected scalar susceptibility. It is related to the two-point 

correlation function according to the formula [^] 

_2 roo ,00 rmrmW-X^RciW) 
Xd V Jo Jo (A 2 + m2)(A 2 + m|)(A' 2 + mf)(A' 2 + m|) { 1 

where -R C (A, A') is the connected two-point correlation function, 

i? c (A,A') = (p(A)p(A'))-(p(A))(p(A')). (14) 



14 



N 
N 
N 



128 
256 
512 



o 

0.001 



0.01 



0.1 



m 



Fig. 7. The valence quark mass dependence of the chiral condensate for 
Nf = and a total number of instantons as indicated in the legend of the 
figure. 

In the derivation of fll3|) we have used that that R c satisfies the sum-rule / °° dXR c (X, A') = 
0. In the bulk of the spectrum, a number variance given by £2(71) — X n is reproduced by 
a two-point correlation function R C (X, A') = x(5(A — A')p(A) — 2p(X)p(X')/N) (For an ex- 
tensive discussion of the relation between a linear dependence of the number variance and 
the sum rule for the two-point function we refer to This results in the disconnected 

susceptibility 

Xd = x(-~ <9m£ - 2£ 2 ) . (15) 
Al/m behavior for the susceptibility was observed in unquenched lattice QCD simulations 



54}|. In Fig. 8 we show the disconnected scalar susceptibility mxd/£(m) versus m. We 



observe a 'plateau' that would indicate al/m dependence for iV = 128, but it disappears 
for larger volumes. Apparently, finite size effects prevent us from making any definite 
conclusions. 

Let us finally make a naive estimate of the disconnected susceptibility. The two-point 
correlation function T 2 (A, A'), which is related to R C (X, A') by T 2 (A, A') = S(X — X')p(X) — 
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Fig. 8. The disconnected chiral susceptibility for Nf 
512 instantons. 



both for 128 and 



R C (X, A'), can be expressed in terms of the two-level cluster function as 

T 2 (A,A')«p(A)p(A')y 2 (r), (16) 
where r = f£ p(s)ds. According to || Y 2 {r) is related to the number variance as 



Yo(r = n) 



(17) 



This results in T 2 (A,A') « i9 A 9 A /S 2 (r). Combining ([16]) and (0) with (0) we find that 
a finite value of Xd, which is expected according to quenched chiral perturbation theory 
49] and lattice QCD simulations |j54| , requires that 

S 2 (A^x) 



lim 



N 



18) 



is finite. In other words, a nonzero disconnected chiral susceptibility requires a number 
variance that shows a linear behavior on macroscopic scales. In chiral random matrix 
theory we find that Xd = in the thermodynamic limit. Because of the large volume 
dependence of the spectral density this argument could not be verified in detail. In future 
work, though, we hope confirm relations of this type. 
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6 Conclusions 



We have identified an energy scale below which the eigenvalue correlations of the QCD 
Dirac operator are given by the chiral random matrix ensembles. In analogy with the 
theory of mesoscopic systems, this scale will be called the Thouless energy For eigenvalues 
near zero the volume dependence of the Thouless energy is weak but its numerical value is 
of the same order as estimates from the pion Compton wavelength. For eigenvalues in the 
bulk of the spectrum the Thouless energy scales roughly with the square root of the volume 
which is in agreement with results from the theory of mesoscopic systems. However, we 
find a proportionality constant that is larger than given by our simple estimates. 

For energy scales beyond the Thouless energy a linear n-dependence of the number 
variance is found which, according to the work by Altshuler and Shklovskii, corresponds 
to a critical system. We have shown that the corresponding Dirac eigenfunctions show 
a multifractal behavior. However, our multicritical exponents are not in agreement with 
relations derived for a critical system. Obviously, our ensemble of instantons is below 
critical. One way to increase the criticality of the system is to increase the number of 
flavors. However, at Nf ^ simulations for the large volumes that are required for 
the extraction of the multifractality index became to costly. Alternatively, one expects 
that pion loops become important at energy scales beyond the Thouless energy, and that 
results from chiral perturbation theory apply to this region. Due to large finite size effects 
we were unable to make a detailed comparison. However, the spectral density near zero 
in the quenched approximation is consistent with a logarithmic dependence predicted by 
quenched chiral perturbation theory. 
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